Effects of Collisional Decoherence on Multipartite Entanglement - 
How would entanglement not be relatively common? 
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Abstract 

We consider the collision model of Ziman et al. and study the robustness of iV-qubit Greenberger- 
Horne-Zeilinger (GHZ), W, and linear cluster states. Our results show that iV-qubit entanglement 
of GHZ states would be extremely fragile under collisional decoherence, and that of W states could 
be more robust than of linear cluster states. We indicate that the collision model of Ziman et al. 
could provide a physical mechanism to some known results in this area of investigations. More 
importantly, we show that it could give a clue as to how iV-partite distillable entanglement would 
be relatively rare in our macroscopic classical world. 



Entanglement lies at the heart of quantum information theory [lj]. It has recently been 
studied extensively in the context of quantum computation, which would require the creation 
and maintenance of highly complex entangled states of TV particles, Ai, • • • , A^, where N 
can be arbitrarily large. This is a difficult task because unavoidable random interactions of 
the composite system A\ • • • A^ with its environment creates entanglement between them, 
which at the same time degrades the entanglement within the system itself. Active interven- 
tion in the form of distillation protocols and error correcting codes is capable of maintaining 
the iV-partite entanglement in the presence of decoherence. However, the efficiency of such 
methods expectedly depends on the a priori robustness of the entangled state in question, 
i.e., the ability of an entangled state to remain entangled under (local) decoherence. Local 
here means that individual particle of the system interacts independently with the envi- 
ronment. Several natural questions of practical importance are thus: What is the robustness 
of multipartite entanglement under certain decoherence model? How the robustness of the 
entanglement changes with the size N of the system? How the robustness of one multipartite 
entangled state compares with other multipartite entangled states? 

In order to gain some insight into the above questions, Simon and Kempe were the 
first to analyze the entanglement properties of iV-qubit Greenberger-Horne-Zeilinger (GHZ) 
states |3[ under the action of (homogeneous) local depolarizing channels. They studied how 
much local depolarization is possible such that the states are still entangled. They made 
comparisons with the W states Q, albeit for W = 3,4 Late, Baadyopadhyay aad Lidar 
p] considered the robustness of iV-qubit GHZ states under local generalized depolarizing 
channels. Carvalho, Mintert and Buchleitner extended the above studies of the iV-qubit 
GHZ and W states to situations involving the presence of local amplitude damping and 
dephasing channels. In 0,0], Dur, Briegel and Hein studied the entanglement properties of 
the rather general class of graph states, which includes GHZ and cluster states [9j], under 
local amplitude damping and generalized deoplarization. 

The structure of iV-qubit entanglement is considerably more complex than bipartite en- 
tanglement. In particular, there is no such thing (yet) as a general practical criteria for 



multipartite entanglement (compare [K L ll|]). So, to analyze the entanglement of the above 
multiqubit states, Simon and Kempe [2| considered the entanglement for bipartite cuts, 
where the N qubits are partitioned into two (arbitrarily chosen) subsystems, with some n 
qubits constitute one subsystem V\ and all the other N — n qubits the other subsystem 
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Vi- For each such cut they employed the partial transposition criterion [h]]. As long as 
the partial transpose of a state has negative eigenvalues, the state is definitely entangled. It 
turns out that a necessary condition for iV-qubit distillability is that the partial transposi- 



tions with respect to any group of qubits are nonpositive (see |l2l Il3||). Consider N distinct 
parties each holding a qubit belonging to the iV-qubit state Pai-a n - It is iV-party distillable 
entangled if there exists a local protocol (i.e., the parties act individually on their qubits 
and are allowed to communicate classically) such that one can obtain from a sufficiently 
large number of copies of pa 1 -a n some true iV-qubit entangled pure state. This distillable 
entanglement criterion is being used in 0,0]. It is of obvious practical relevance and will be 
the criterion employed in this paper. We remark that Carvalho, et a/.[6( too used a different 
measure of multipartite entanglement - a specific generalization of concurrence 14] of pure 
bipartite states, that is sensitive to real multipartite correlations. 

Carvalho, et al. .6:] showed that only for depolarizing channels do the GHZ or W states 
turn separable after a finite time. And, this occurs with the concurrence of the decaying 
GHZ state always larger than that of the decaying W state. In contrast, the amplitude 
damping as well as the dephasing channels induce separability only in the limit of infinite 
time. And, the concurrence of the decaying W state gets larger than that of the decaying 
GHZ state after a short time. In terms of entanglement decay rates, they found that the 
GHZ states become separable at rates which increase linearly with N (for iV > 4). For the W 
states, only the depolarizing channels give rise to an almost linear increase of entanglement 
decay rates with N, slightly faster than for the GHZ states. The decay of the concurrence is 
independent of N for amplitude damping and dephasing channels. They concluded that the 
multipartite quantum correlations of W states outperform those of GHZ states in terms of 
their robustness. Their results for the GHZ states agree with those in [?], 0] , which concluded 
via their analysis on linear cluster and more general graph states that true multipartite 
entanglement in macroscopic objects can be more stable and might be more common than 
previously thought 1151 1 . 

nn 

Recently, Ziman et al. [la, 113] showed that all decoherence maps of qubits, to be defined 
below, can be modelled as sequences of collisions between the system qubits under consid- 
eration and qubits from their environment [3]. In this paper, we study the robustness of 
multiqubit GHZ, W and linear cluster states under collisional decoherence, which we briefly 
describe next. We find that the iV-qubit distillable entanglement associated with these states 
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may not be robust, with the GHZ states the most fragile, and the linear cluster states less 
robust than W states. 

Consider a system of N qubits, A ly ■ ■ ■ , A N , initially decoupled from an environment, 
i.e., the initial system-environment state is given by pa\—a n ® S era) . Following [16], we 
model the environment set of nili N[ qubits: E tj , with 1 < i < N and 1 < j < N[. 
These environment qubits are initially in a factorized state: E env = (8)^1 1 Ce„ j with 
^E i:j — £ f° r an Eij- They do not interact between themselves, and each E^ undergoes a 
bipartite collision with each system qubit Ai only once. Under such conditions the evolution 
of each Ai is described by the sequence of maps £ Ai , £ A .. In particular, the state 

of Ai after the fcj-th collision is given by p A *' = S^.IS^' 1 ^ ■ ■ S A .[p Ai ] ■ ■ ■]], where p Ai = 
tr Al ... Ai _ lAi+1 ... AN (p Al ... AN ), S J A .[p Ai \ = tr Ei .[U A . Eij (p Ai ® £ Erj )U AiEij ]. Here, the bipartite 
unitary operator that describes each collision, 

U AiEij = |0)a,(0| ® vl + |l) Ai (l| ® v x Ei , (1) 

It belongs to the set of all controlled-f/ transformations, where the system qubit Ai plays 
the role of the control and the environment qubit E^ is a target. With the unitary operators 
on Eij defined by 

where \tp) and \<p) are arbitrary normalized kets with (ip^-fy) = (^{(p) = 0, it follows that 
the state of A{ after the first collision is given by p^ A . or 

EilPA*) = PaA + (x Eil hP%sX + (x Ei1 )spIs a > + PaA (3) 

where 




S + = \0)(l\= I ° * |,5- = |1><0|= I ; : |, (4) 

such that p Ai = p%P% + P%S% + P Ai S Ai + P A \P Ai - Here, the unitary operator X E .. = 
V E \yl., (X Ei .)s= (X Eij )l and 

(Xfiyjf = tr[£X £i .] ee Ay exp(/<^), (5) 
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with < Ay < 1 and < </>y < 2tt. Whenever the collision is nontrivial, i.e., V^_. ^ 
e ~i<Pij yo ; anc i the environment qubit state £ is not an eigenstate of , we have < \j < 1 
and Eq.(3) describes a quantum decoherence channel with decoherence basis {|0), |1)} 
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17(. From here on, we assume these conditions are fulfilled. Specifically, after the /cj-th 



collision, we have P^ 1 — > P^. , 

S±-^ lt eM±I& t )SX, (6) 

where 

hi hi 

7i = Y[X iv = Ys&r ( 7 ) 

i=i j=i 

The most nontrivial j-th collision occurs when Ay = 0, and it results in 7^ = 0. If A\ does 
not collide with any environment qubit, then 7, = 1 and $j = 0. We remark that the above 
model is intrinsically discrete and has not been considered for analysis of robustness of N- 
partite entanglement. However, Ziman et al. have shown that provided < Ay < 1 one can 
write down the master equation for a dephasing channel by performing some continuous- 



time approximation 
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17l |. For nonzero Ay's, the model thus could provide a microscopic 
description of the mechanism for a dephasing channel. More significantly, we show that 
the model indeed yields new results quite different from previous investigations, especially 
when Ay = 0. It is beyond the scope of this paper to discuss exactly how probable it is for 
Xij = to occur. However, only looking at Eq.(5) superficially, this may just happen often: 
for instance, when £ = a°/2 and X = sin 9 cos (pa 1 + sin 9 sin (pa 2 + cos^a 3 , where a° is the 
2x2 identity matrix, cr l 's are the Pauli matrices, < 9 < it, and < <p < 2tt. 

We start out by investigating the entanglement properties of iV-qubit GHZ states, 
\^ghz)a 1 -a n = (|0 ■ ■ ■ 0)a 1 -a n + |1 ■ • • 1)ai-Ajv)/\/2 under this collisional decoherence. 



\^ghz)a!- a n has the density matrix 

N N N N 

lv 2 ' 



PTX = I [ F% + 31 + (g) Si + Pi ) . (8) 

i=l i=l i=l i=l / 



Applying the decoherence channel, Eq.(3), to every system qubit yields 

| N N N N "> 

i.v _ 2 



pm z a n = I [ (g) f% + r <g> si + r ® + <g> Pi ) , (9) 

i=l i=l i=l i=l J 



with T = (iliLiTij ex P [iJliLi^i)- Eq.(9) follows from Eq.(6). Take partial transposition 
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of the first n qubits. Since P4. 1 — > P%. and S^. — > fij., we have 

N 

(10) 

Let X = (p'g HZ AN ) Tvi -*Im.-a n = {i~x) ®f =1 P%-x®U Pl®P\ N x Pi ® 

+ (I - x) ®ti ^1, + r OLi fi£ <8>iL n+1 ^ + T* <g)? =1 fij ®iL„ + i fi^ • det * = then 
yields x 2N ~ 4 (l — x) (x 2 — ||r| 2 ) = 0. It is not difficult to deduce that, with respect to all 
bipartite cuts, {.p'a^ Z a n ) Tvi nas the same negative eigenvalue given by 

1 AT I N k t 

j^5j=~n7<=~nnv ( n ) 

z t=i z 1=1 j=i 

We note that a computable measure of entanglement associated with a bipartite cut V1WV2, 
the negativity Q of pg?.^ (or pgjf£) is given by lA^pg^JI- If Ay = 1 for all i and 
j, then \N{p'g HZ AN }\ = 1/2- As long as < Ay < 1, |Af [p^l J | is nonzero, i.e., pg^ w 
remains iV-partite entangled. Obviously N[p'a^. Z a N ) depends on both N and hi. When 
hi = K and Ay = A for all % and j, we have 

X[pl HZ A N ] = ~\^ KN - (12) 

n 

Eq.(12) agrees with the results of analysis for the dephasing channel in jjy . For fixed A and 
N, \N[p'a^- Z a J I § oes exponentially to zero as K — > 00, and for fixed K and A, In |A/"[p^f.^ ] | 
varies linearly with N. However, we emphasize here that all it takes is merely one Ay = to 
reduce ^[p'a^An] t° 0, and in this case, for p'a^ Z a n to become separable. In this sense, we 
conclude that the iV-qubit entanglement associated with a GHZ state is extremely fragile 
under collisional decoherence. 

Next, we turn our attention to the case of iV-qubit W states, \^w)a 1 ---a n = 
(|0 ■ ■ • 01) Ai - a n + |0 ■ • • 010) a 1 -a n + - ■ • + 1 10 • ■ ■ 0)a 1 -a n )/VN, which has the density matrix 
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N 



PAi-A n E + ^2 ^mrn' + ^mm' > (13) 

iV \m=l l<m<m'<N ) 

with H m ee Pi ® P\ m ® ®«= m+ i n and T^ m , = ®T=x P° Ai ® ^1 ® ^ ® 

fi^ , ® ^ilm'+i-^i- There are altogether iV(iV — 1) terms in the second summation. To 
any term with a pair S\ m and ,, there is a corresponding term with the pair S^ m and 
S~a , ■ Applying the decoherence channel, Eq.(3), to every qubit yields 

AT 



JW 1 



N 



PAx-A 



mm' mm 

m=l l<m<m'<N 



(14) 
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where T mm / = jmlm' exp[J(<l ) m — $ m ')]. As above, taking partial transposition of the first 
n qubits, we have 

N 



(JW \T Vl _ J_ 

\Pa v -a n ) - N 



71=1 l<m<m'<n 
l<m<n<m'<A r 



(15) 

n+l<m<m'<A r 

where = (g)^ 1 1% ® Sj m ® <g)™Ui P£ ® S£„ ® 1%. Similarly, we consider 

X = (p'a V 1 ---a n ) Tvi ~ x Ia 1 ~ a n and det x — 0. However, in contrast to Eq.(ll), depending on 
the bipartite cuts, {p'Y 1 -A N ) Tvi has negative eigenvalue given by 



WZ..aJ = ~ 



, . E 7f E 7?|. (16) 



And, the negativity associated with a bipartite cut Vi\\V% is given by lA^^...^ ]|. If Xij = 1 



for all z and j, then |A/"[p^... Ajv ] | = yn(N - n)/N. Eq.(16) agrees with Eq.(ll) for N = 2. 
As above, remains entangled so long as < A^ < 1. However, M[p'Y r .. A ] clearly 

has a very different dependence on both iV and fcj compared to A/"[p^^ ]. When ki = K 
and A,j = A for all z and j, we have 

A^K-aJ = —y/niN-n)* 1 *, (17) 

which reduces to —X 2K /2 independent of N if n = N/2. On the other hand, when only i — 1 
and fci = i^, we have either 

1 



or 

JW 



Wt..A N ] = -jjJn(N -n-l + X 2K ), (19) 

depending on if Ai G V\ or G TV These remain nonzero even in the limit of infinite 
K. Hence, the "weakest links" come from Eq.(17_). Again, Eq.(17) is consistent with the 
results of analysis for the dephasing channel in p. For fixed A and N, \N\p'Y v ..A n ]\ § oes 
exponentially to zero as K — > oo. \J\f[p')^...A J I decreases as the number of qubits N increases. 
But, this decrease is solely due to the property of the original pure W state. Therefore, 



modulo this effect, the entanglement decay rate is independent of N. Here, we emphasize 
that in contrast to Eq.(ll), Eq.(16) implies that when A y - = or 7$ = for one A{, the state 
of the remaining N — 1 qubits A 1; • • • , • • • , Ajv is still (AT — l)-partite distillable 

entangled, even though p'Y v ..a n ls 110 l° n g er Af-partite distillable entangled. Furthermore, it 
has to take ji = for all to reduce the Af-qubit state to a completely separable one. We 
may thus conclude that the multipartite correlations of W states outperform those of GHZ 
states in terms of their robustness under collisional decoherence. 

Lastly, we consider a linear cluster state of A" qubits . The state can be written in the 
form 

HT N 

\*c)m~a» = <8>(l°U4 i+1 + |1>*) (20) 



i=i 



with the convention &\ N+i = 1. For N = 2, applying the docherence channel, Eq.(3), to 
both qubits, we obtain the resulting state p'a x a 2 with negativity 

£[Pma 2 \ =max{^ 12 ,0}, (21) 

where r/12 = (7172 + 7i + 72 — l)/4. When A y - = 1 for all i and j, it gives 1/2, which 
agrees with both Eqs.(ll) and (16). However, for < Ay < 1, it obviously yields different 
result, This is despite the tact that \* cUlA , is up to a local unitary operation a nraxnnally 
entangled Bell state Consequently, supposing 71 = 72 = 7, the bipartite state p'a 1 a 2 
becomes separable at 7^ = — 1 + \pl. In this sense, the bipartite linear cluster state is less 
robust compared to the corresponding GHZ or W states. This fact that the robustness of 
a state can be affected by local operations has been discussed in |2j. Similarly, for A^ = 3 
we have Pa 1 a 2 a 3 w hh negativity associated with the bipartite cut Ai e V\ and A 2 ,A 3 e V2 
given by 

^[PA 1 (A 2 A 3 )] = max{r /l2 ,0}. (22) 
By symmetry, we have with 7/23 = (7273 + 72 + 73 — l)/4, 

£[P(^ 2 )J= max te3,0}. (23) 

But for the bipartite cut Ax, A 3 G Vx and A 2 we have instead 

^[p'( C A 1 )A 2 (A 3 )} = max{r/ 12 , ^23, V123, 0}, (24) 

where 77123 = [(1 + 71)72(1 + 73) - (1 - 7i)(l ~ 73)]/8. Eqs.(24) - (26) give 1/2, which agrees 
with Eq.(ll) when \j = 1 for all i and j. This is again due to the fact that a tripartite 
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linear cluster state is up to local unitary operation equivalent to a tripartite GHZ state |9j . 
For < Xij < 1, and assuming 7l = 72 = 73 = 7, £\pf AlM)A ^ = £[Pa i{ a 2 a 3 )} = at 7^ but 
^[p'{a 1 )a 2 (a 3 )\ — at 7^?] t ~ 0.295598. We again conclude that the tripartite linear cluster 
state is less robust compared to the corresponding GHZ state. The above cases exhibit two 
important properties generic to all iV-partite linear cluster states. First, among the negative 

eigenvalues of the partial transposed states across all possible bipartite cuts V\ and V2 are 

(3) (2) 

Vijj Vijk, etc. Second, although %rlt < IcAt it is nonzero, and the weakest links come from the 
%'s. Unlike the W states, \$ c) a x -a n can thus fail to be iV-partite distillable entangled for 
a range of nonzero values of 7^. Therefore, the linear cluster states are less robust compared 
to the W states, even though in this case it also has to take 7« = for all A, to reduce an 
iV-qubit linear cluster state to a completely separable one. 

In conclusion, we have shown that under collisional decoherence the multipartite en- 
tanglement associated with the GHZ, W and linear cluster states may not be robust. In 
particular, an iV-partite state can easily fail to remain iV-partite distillable entangled. Our 
results may seem contrary to some recent investigations. However, it is not surprising, as 
this is the first time an inherently discrete model of decoherence has been considered for 
analysis of this kind. It shows that the debate over if multipartite entanglement should be 
relatively common is far from over. We end with the following remark: should there be 
some sort of mechanism that causes Ay = |tr[^Xs y ]| = to occur very frequently, it could 
provide an answer to why multipartite entangled states are not prevalent in our macrscopic 
classical world. 
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